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^ , The energy dependence of the radii size of jets are studied in detail by the HBT cor- 

(-H 1 relation method using Monte Carlo Simulation generator Jetset7.4 to produce 40,000,000 

■ ■ events of e+e" collisions at = 30, 50, 70, 91.2, 110, 130, 150 and 170 GeV. The radii 

of jets are measured using the HBT correlation method with the indistinguishability of 
identical final state pions. It is found that the average radii of quark-jets and gluon-jcts 
are independent of the cm. energy of e+e~ collisions. The average radius of quark-jets 
are obviously larger than that of gluon-jets. The invariable average radii of quark-jets 
and gluon-jets in e+e~ collisions are obtained at the end of parton cvolvcmcnt. 
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CO ' 1. Introduction 

It was well-known that Hanbury-Brown and Twis had brought forward HBT cor- 
relation in the process of measuring the angular radii of the emitting sources in 
Later, the HBT correlation was widely used in subatom studies - ^l*^. 
The HBT correlation method has been an important way to measure the size of the 
emitting source in high energy collisions. Due to "color confinement", we cannot 
observe free quarks and gluons and cannot yet measure the size of them directly. 
However, the HBT correlation method offers a viable indirectly method, and apply- 
ing this method into the high energy collisions we can obtain some characteristics 
of strong interaction for quarks and gluons. 

Historically, the discovery in 1975 of a two-jet structure''' in e+c^ collisions at 
center of mass (cm.) energy > 6 GeV had been taken as an experimental confir- 
mation of the parton modepl, and the observation in 1979 of a third jet in e'^e" 
collisions at 17 — 30 G eV had been recognized as the first experimental evidence of 
the gluon I 9 | f0 | f f I S I ^j^e early 199 Qs, the production of jets in hadron-hadron 
collisions was widely studied I J^3 | 14 | 15 | Yia,d been considered as an efficient way 

to obtain the strong coupling constant as How to distinguish jets and the study 

1 7 1 S 

on jets are also very important, in relative high energy ion collisions . Based 
on this idea we can get information about quarks and strong interactions from the 
study of jets by using the HBT correlation method. 



1 
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In e"'"e~ collisionil^, firstly, the e+e" pair is annihilated into a virtual 7*/Z° 
resonance. The virtual 7*/^", in turn, decays into a qq pair. Then the initial qq 
may radiate other gluon and qq pairs, giving rise to a cascade process. This stage 
is responsible for the formation of hadronic jets. Further, the unstable hadrons 
decay into experimentally observable particles (mostly pions). It has been found 
that the majority of e^e~ collision events have a 2-jet structure. If an initial quark 
or anti-quark emits a hard gluon with sufficiently large transverse momentum, a 
3-jet structure can be formed. Thus, the source of a single jet is from a single initial 
quark (or anti-quark) or gluon. 

Although the quark and gluon, before being observed, have been fragmented into 
the final state hadrons, the final state particles inside the jets still carry a lot of in- 
formation about the parent quark and gluon. The quark and gluon are two different 
types of particle. For example, the quark is a fermion with colour charge equal to 
4/3, while the gluon is a boson, carrying colour charge 3. These differences will cer- 
tainly influence their fragmentation, resulting in different properties of quark-jets 
and gluon-jets.Some characteristics of quarks (anti-quarks) or gluons is reflected by 
the geometrical characteristics of jets. So, the study of the geometrical characteris- 
tics of the jets is helpful in the understanding of the perturbative/nonperturbative 
properties of QCD. 

I on I 

In the vefy^, the geometrical characters of quark-jets and gluon-jet have been 
studied with the HBT correlation method using MC generator producing quark- 
jets and gluon-jets in 3-jet events of e+e~ collisions at ^/s = 91.2GeV. However, do 
the size of quark-jets and gluon-jets depend on the cm. energy of e"''e~ collisions 
producing these jets? Are the size of quark-jets measuring in 3- and 2-jets events 
of e+e^ collisions the same? Our work will focus on these questions. 

The paper is organized as follows: In Sec. II, we briefly introduce the method of 
identiflcation jets and the HBT correlation function. In Sec. HI, the average radius 
of quark-jets and gluon-jets in 2-jet Events are calculated. In Sec. IV, the average 
radius of quark-jets and gluon-jets in 3-jet Events are calculated. A short summary 
is the content of Sec. V. 



2. The method of identification jets and HBT Correlation 
Function 

In our work, the data of e+e^ collision events are produced by Monte Carlo Simu- 
lation generator Jetset7.4. The 2-jet events and the 3-jet events are selected using 
the Durham jet algorithirJ^. In these methods, there is a cutting parameter j/cut, 
which, in the case of the Durham algorithm, is related to the relative transverse 
momentum kt aJ^ 

h = y/Vcut ■ \fs, (1) 

where is the cm. energy of the collision. From the experimental point of view, 
kt can be taken as the transition scale between the hard and soft processes. Its 
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value depends on the definition of "jet" . 

Tlie sin gle q uark-jet and single gluon-jet are identified from 3-jet events using the 

OOI 

angular rule^ .We assume that the three jets in a 3-jet event come from quark, anti- 
quark and gluon, respectively. Because of energy-momentum conservation, the three 
jets in one event must lie in a plane, which is shown in Fig.l, where Pi{i = 1, 2, 3) 
is the total momentum of all particles in jet-i. The jets are tagged using the angles 
between them: 



P3 jet 3 




jet 2 



Fig. 1. The sketch map of the 3-jet distribution. 



arccos 



PjlPkl+Pj2Pk2+Pj3Pk3 



,VPh+Ph+PhVPL+Pi2+pLj ' (2) 

{i,j,k = 1,2,3; i 7^ j;j ^ k;k^ i), 

where the largest angle, 6*3, faces the gluon jet; the smallest angle, di, faces the jet 
formed by an initial quark without emitting a hard gluon; and the middle one, 02, 
faces the mother-quark-jet. 

According to the requirement of momentum conservation, the three jets should 
be in one plane, and we add the condition: di + 62 + O3 > 359° . 

It should be noticed that the angle opposite to the mother-quark-jet is very 
close to that opposite to the gluon-jet, i.e. 62 ~ ^3, so they are easily confused. 
Therefore, we demand a cut condition: 93 — 62 > Ocut , here 9cut = 10°. This cut 
rejects about 12% eventi^S. 

The HBT correlation, also called the Bose-Einstein correlation, results from the 
indistinguishability of identical final state particles. Most of the final state particles 
produced in e^e^ collisions are tt mesons, so we choose tt mesons (7r+,7r~, vr") as 
the identical particles to study. If P{ki, k2) is defined as the probability of observing 
two identical pions at the same time with momentum fci and fc2, and -P(fci) and 
P{k2) are defined the probability of observing pions with the momentum fci and 
^2, respectively. The correlation function (72(^1,^2) is defined as: 

C2{kuk2) ^ p^^^^p^^^y (3) 

If the equivalent density function of the source is parameterized to Gaussian 
form, we have: 

C2(fci,fc2) = C2(g, fci,fc2) , . 

= 1 + A exp{-i?2q2 „ ^2^2 „ ^2^2 _ ^2^2}^ W 
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where g = fci — fc2 is the four-dimensional momentum difference. 

If only the spatial part of the source is considered and assume that the dis- 
tribution of the source is isotropic, the correlation function can be simplified as: 



C2{ki,k2) = C2{q,ki,k2) 

= 1 + aexp{-i?2g2| 



2^21 (^) 



where ^ = k i — k 2 '^s the three-dimensional momentum difference. According to 
the definition, jets do not possess spherical symmetry, but are axially symmetric 
instead. So the parameter R characterizing the geometrical properties of the jets is 
actually the average radius of jets. 

In this paper, we study the average radius R of the pion source only through the 
spatial distribution function of the source, which is taken as spherically symmetric. 
Then, the information about the average size of the emitting source for the final 
state particles can be obtained. 

The three-dimensional momentum interval region chosen is Q = ^ 2.5 GeV/c, 
and is equally divided into 50 cells. We use Monte Carlo simulation generator 
Jetsct7.4 to produce e+e~ collision events both with and without HBT correlation, 
and then select out suitable events for study. Identical tt mesons are selected from 
the final state particles to make pion pairs after any two tt mesons are grouped 
with each other. The three-dimensional momentum difference of the tt meson pairs 
are calculated. The correlation function (also called correlation coefficient) with 
statistical method is: 

C2{k„k2)=C2{Q)=^, (6) 

where Fc{Q) is the three-dimensional distribution function of the identical particles 
with HBT correlation inside the jets and F{Q) is the three-dimensional distribution 
function of the identical particles without HBT correlation inside the jets. Since the 
correlation among identical particles with large momentum difference is quite weak, 
the distribution here with the HBT correlation should be almost the same as the 
distribution without the HBT correlation. Thus the C2{Q) can be multiplied by a 
coefficient to make the value of it equal to 1. Thus, using Eq(5) to calculate the 
average radius R, the Eq(6) can be expressed as: 

C2(0) = ^ = r,(l + acM-R'Q')), (7) 
where r] is the value of the correlation function C2 (Q) at a large momentum interval. 



3. The Measurement of source radii inside jets of 2-jet events 

We use Monte Carlo Simulation generator Jetset7.4 to produce 40,000,000 events 
of e+e" collisions both with and without HBT correlation, which the cm. energies 
are ^/s = 30, 50, 70, 91.2, 110, 130, 150, 170 GeV,respectively. The final state tt 
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mesons (7r+,7r~, 7r°) as the identical particles are chosen from event samples to 
study. We just select out the final state identical n mesons emitted from vertex at 
origin, because some of them are secondary emitted or even multistage emitted. 
Thus, the calculated result is able to reflect the characteristics of the source of jets 
properly. The 2-jet events are selected using the Durham jet algorithm, and the 
cutting parameters at different cm. energies are selected as Table 1. 



Table 1. The cutting parameters of selection 2-jet events from e^"e colli- 
sion events. 



(GeV) 


30 


50 


70 


91.2 


110 


130 


150 


170 


2/cut 


0.07 


0.06 


0.06 


0.05 


0.078 


0.078 


0.078 


0.078 



The 2-jet event is constituted by the two jets which arc formed by the fragmen- 
tation of the back-to-back qq. And the two jets are called quark jet-1 and quark 
jet-2, respectively. Due to the back-to-back symmetry of the two jets formed by the 
fragmentation of the qq, the distribution patterns of the two jets should be totally 
the same. So, we just need to study one of the two jets. We will choose quark jet-1 
which is referred to as quark-jet. The correlation function is produced both with 
and without HBT correlations when the three-dimensional momentum interval re- 
gion of the identical tt mesons is chosen as Q = ^ 2.5 GeV. According to equation 
(6) we calculate the value of the correlation function C2(Q) of tt (7r+,7r~ and tt") 
mesons inside quark-jets from 2-jet events for all the 8 cm. energies. Then, the 
average radii size of emitting source of pion meson inside single jet can be obtained 
through fitting the correlation functions C2{Q) using Eq.(7) for 7r"'",7r^ and tt" , 
as shown in Table 2, respectively. As an example, the results for 3 cm. energies are 
shown in Fig. 2 

It is clear to see from Fig. 2 that the distributions of correlation functions of 
the TT meson inside quark-jets at different cm. energies are similar. Especially, the 
distributions of particle both for 7r+ and 7r~ mesons at the same energies are in 
superposition. However, the distributions of identical particles for tt" mesons is 
different to 7r+ and 7r~ mesons owing to the electromagnetic interaction among 7r+ 
and 7r~ mesons in the process of hadronization. The mean radii R of jets for all 
various cm. energies are listed in Table 2 at the last line. 

For the convenience of comparison, we draw the source radii of the three kinds 
of TT mesons inside quark-jets for all energies in Fig. 3. 

It can be seen from Table 2 or Fig. 3 that the values of the radii R of jets for one 
meson from different cm. energies are nearly the same within the error range; the 
radii R of jets for the same cm. energy both 7r+ and 7r~ mesons are approximately 
the same within the error range, i.e. their means Rq,Tz+ = 0.693 ±0.003 fm, Rq^T^- = 
0.690 ± 0.004 fm; But there are some difference between the 7r+ or tt" mesons and 
7r° mesons for the Rq^T^o = 0.6000 ± 0.006 fm. This is due to the electromagnetic 
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Fig. 2. The distributions of correlation functions of quark-jet from 2-jet events at different cm. 
energies. The values of cm. energies are (a) 50 GeV, (b) 91.2GeV, (c) 150 GeV, respectively. The 
data for the e+e~ collisions are produced by MC Jetset7.4 generator. 



Table 2. The radii Rg of quark-jets from 2-jet events at 
different cm. energies, measured by HBT using jr+jTr^ 
and 7r~ mesons, respectively. 



V5 (GeV) 








30 


0.704±0.006 


0.553±0.007 


0.708±0.006 


50 


0.662±0.004 


0.580±0.006 


0.686±0.004 


70 


0.703±0.004 


0.577±0.005 


0.688±0.004 


91.2 


0.672±0.003 


0.607±0.006 


0.684±0.004 


110 


0.697±0.003 


0.593±0.004 


0.678±0.004 


130 


0.710±0.004 


0.616±0.005 


0.679±0.003 


150 


0.693±0.005 


0.624±0.007 


0.685±0.005 


170 


0.705±0.004 


0.650±0.008 


0.715±0.004 


-Raver 


0.693±0.004 


0.600±0.006 


0.690±0.004 



1 
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Fig. 3. The comparison of the values of the radii R of quark-jets from 2-jet events for different 
types of TT mesons (7r+,7r~, n^) at different cm. energies. 



interaction among 7r+ and n mesons in the process of hadronization. So, the results 
for TT^ mesons are more authentic. 



20, 2013 2:29 WSPC/INSTRUCTION FILE Chen-201302 



7 

4. The Measurement of source radii inside jets of 3-jet events 

In the same way, we use Monte Carlo Simulation generator Jetset7.4 to produce 
40,000,000 events of e+e~ collisions both with and without HBT correlation for 8 
cm. energies. And we choose the final state tt mesons {tt^,tt~, 7r°) emitted from 
the vertex at origin as the identical particles to study. The 3-jet events are selected 
using the Durham jet algorithm, which the cutting parameters for 8 energies are 
listed in Table 3. 



Table 3. The cutting parameters of selection 3-jet events from e+e collision events. 



(GeV) 


30 


50 


70 


91.2 


110 


130 


150 


170 


Vcut 


0.008 


0.005 


0.004 


0.002 


0.0015 


0.001 


0.0008 


0.0005 



After the quark-jet, mother quark-jet and gluon-jet arc identified from 3-jet 
events using the angular rule, the correlation function is produced both for the 
case with and without HBT correlation when the three-dimensional momentum 
interval region of the identical tt mesons is chosen as Q = ^ 2.5 GeV/c. We 
calculate the value of C2{Q) according to formula (6), then the average radius size 
of emitting source of pion inside single jet can be obtained through fitting the 
correlation functions C2(<3) using Eq.(7) with 7r+,7r~ and 7r° meson, respectively. 
The results are listed in Table 4. As an example, the results for 3 cm. energies are 
shown in Fig. 4. 

It is easy to come to the conclusion from Fig. 4 that the three type of tt meson 
correlation functions for quark-jcts, mother quark-jets and gluon-jcts are all sim- 
ilar. And the distributions of different identical particles, for 7r+ and 7r~ mesons, 
formed at the same energy are nearly in superposition. However, the distributions 
of identical particles for 7r° mesons are different to 7r+ and tt" mesons. 

For the convenient of comparison, we draw the average source radii of the three 
kinds of jets at diflferent cm. energies for the three types of tt mesons in Fig. 5. 

It can be seen from Table 4 and Fig. 5 that: the values of radii R of quark-jets or 
mother quark-jets and or gluon-jets at different cm. energies are all nearly the same 
within the error range, respectively; the average radii of quark-jets are obviously 
larger than that of gluon-jets, and the average radii of quark-jets are also larger 
than that of mother quark-jets. The mean radius of quark-jets for 7r+ mesons is 
0.71 ± 0.02 fm, that of mother quark-jets is 0.67 ± 0.01 fm, and that of gluon-jcts is 
0.59 ± 0.01 fm. The mean radius of quark-jets for 7r° mesons is 0.61 ± 0.02 fm, that 
of mother quark-jets is 0.55 ± 0.02 fm, and that of gluon-jets is 0.43 ± 0.01 fm. The 
mean radius of quark-jets for tt" mesons is 0.74 ±0.01 fm, that of mother quark-jets 
is 0.68 ± 0.02 fm, and that of gluon-jets is 0.60 ± 0.01 fm,respectively. 
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Fig. 4. The distributions of tt meson correlation functions of the three kinds of jets from 3-jet 
events at different cm. energies. The first row is for the case of quark-jets; the second row is for 
the case of mother quark-jets; the third hue row is for the case of gluoii-jets. The cm. energies of 
(a), (d) and (g) are 50 GeV; (b), (e) and (h) are 91.2 GeV; (c), (f) and (i) are 150 GeV,respectively. 
The data for the e+e~ collisions are produced by MC Jetset7.4 generator. 
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Fig. 5. The comparison of the values of the radii R of the three kinds of jets from 3-jet events at 
different cm. energies for different tt mesons, (a) 7r+, (b) tt", (c) tt". 
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Table 4. The values of the source radii R of the three types of jets 
at different cm. energies calculated with the three kinds of tt mesons 
(7r+, 7r~, tt"). 



TT mGSons 


\/s (GeV^ 


''quark V'-**-'-/ 


-'tm — quark 


''gluon V,-'"'-**^ 




30 


0.60±0.03 


0.63±0.02 


0.69±0.02 




50 


0.72±0.02 


0.64±0.02 


0.66±0.01 




70 


0.73±0.01 


0.70±0.01 


0.54±0.01 


7r+ 


91.2 


0.71±0.01 


0.68±0.01 


0.64±0.01 




110 


0.74±0.01 


0.64±0.01 


0.55±0.01 




130 


0.74±0.U1 


O.DoiO.Ul 


O.oDiO.Ul 




150 


0.73±0.01 


0.68±0.01 


0.56±0.01 




170 


0.75±0.01 


0.70±0.02 


0.55±0.01 




^avcr 


0.71±0.02 


0.67±0.01 


0.59±0.01 




30 


0.44±0.02 


0.41±0.02 


0.48±0.02 




50 


O.eOitO.Ol 


0.54±0.02 


0.46±0.01 




70 


o.eiito.oi 


0.53±0.01 


0.42±0.01 




91.2 


0.67±0.02 


0.53±0.01 


0.44±0.01 




110 


0.64±0.02 


0.60±0.01 


0.42±0.01 




130 


0.60±0.01 


0.57±0.02 


0.41±0.01 




150 


0.67±0.02 


0.61±0.02 


0.41±0.01 




170 


0.66±0.02 


0.59±0.02 


0.44±0.01 




^avcr 


0.61±0.02 


0.55±0.02 


0.43±0.01 




30 


0.76±0.02 


0.65±0.03 


0.67d=0.02 




50 


0.75±0.02 


0.66±0.01 


0.62±0.01 




70 


0.72±0.01 


0.70±0.02 


0.60±0.01 


n~ 


91.2 


0.75±0.01 


0.65±0.02 


0.61±0.01 




110 


0.73±0.01 


0.67±0.01 


0.55±0.01 




130 


0.75±0.01 


0.69±0.01 


0.57±0.01 




150 


0.74±0.01 


0.67±0.01 


0.60±0.01 




170 


0.72±0.01 


0.72±0.02 


0.59±0.01 




^aver 


0.74±0.01 


0.68±0.02 


o.eoito.oi 



5. Conclusion and discussion 

We use Monte Carlo Simulation generator Jetset7.4 to produce the data of e+e~ 
collision events, which the cm. energies are ^ = 30, 50, 70, 91.2, 110, 130, 150, 170 
GeV. The 2-jet events and 3-jet events are selected using the Durham jet algorithm. 
The geometrical characters of quark-jets and gluon-jets are studied in detail using 
the HBT correlation method. The conclusions are as follows: 

• The radii of quark-jets or gluon-jets measured at different cm. energies are ap- 
proximately the same within the error range, which shows that the radii of quark- 
jets and gluon-jets reflect some intrinsic properties of quarks and gluons. 

• The values of the mean radii of quark-jets, mother quark-jets and gluon-jets 
are measured by calculating the final state identical tt mesons using the HBT 
correlation method are obtained, shown in Table 5. 

• The mean radii of quark-jets measured is obviously larger than that of gluon, 
which indicates that the size of quark is larger than that of gluon. However, the 
mean radii of mother quark-jets measured is less than that of quark-jets. This 
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Tabic 5. The means radii R of jets measured using the HBT correlation 
method in 2-jct events and 3-jet events from e+e" collision events. 



TT mesons 


2-jet events 




3-jet events 






^quark (fm) 




^im-quark(fm) 


-Rgluon (fm) 


7r+ 


0.693±0.004 


0.71±0.02 


0.67±0.01 


0.59±0.01 




0.600±0.006 


0.61±0.02 


0.55±0.02 


0.43±0.01 


7r~ 


0.690±0.004 


0.74±0.01 


0.68±0.02 


0.60±0.01 



may be due to the mixture o£ a small amount of gluon-jets and mother-quark- 
jets in the process of measurement which makes the radii of mother quark-jets 
measured smaller. 

• The results for tt^ mesons are more authentic than for 7r+ and tt~ mesons for 
there are no electromagnetic interactions among tt*^ mesons in the process of 
hadronization. 
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